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CN . Abstract 

C\| 0{as) QCD corrections to the inclusive B Xse~^e~ decay are investigated within the two - Higgs 

^ ■ doublet extension of the standard model (2HDM). The analysis is performed in the so - called off 

V/-^ ', -resonance region; the dependence of the obtained results on the choice of the renormalization scale 

^ is examined in details. It is shown that 0{as) corrections can suppress the B — > Xse~^e~ decay 

^ _ width up to 1.5 -7- 3 times (depending on the choice of the dilepton invariant mass s and the low - 

^ ' energy scale fi). As a result, in the experimentally allowed range of the parameter space, the relations 

^ between the B XsC^e" branching ratio and the new physics parameters are strongly affected. It 

is found also that though the renormalization scale dependence of the B XgC'^e^ branching is 

I \ significantly reduced, higher order effects in the perturbation theory can still be nonnegligible. 
Q-i. 
D . 
^ ■ 

> : 
•1— I , 

X- 

^3: 



1. Rare B - meson decays can serve as an important source of information on new physics beyond 
the standard modeL The first experimental evidence for these decays has been observed by CLEO 
[0 for the exclusive B K*'j channel. Later on the branching ratio of the inclusive B Xgj decay 
has been measured by CLEO, ALEPH and Belle collaborations [0-0- The weighted average for this 
decay branching is 

^exp^^ ^ X,7) = (2.96 ± 0.35) x 10"^ (1) 

Recently the evidence for the rare exclusive channel B — > Kfi'^fi" has been also observed 0. 

The experimental result for the B Xg'j branching is in a good agreement with the SM pre- 
dictions (see and references therein). The new physics contribution to this decay width must be 
small enough to avoid the contradiction with the experiment, hence studying B —>■ Xg'j decay one 
can get some constraints on the new physics parameters. 

The another popular inclusive rare B decay mode, B Xse~^e~, has not been observed yet 
(only upper bound on its branching ratio exists 0), however it is expected to be measured dur- 
ing the forthcoming experiments at the B - factories. Then, analogously to i? ^ -^s7, the study 
of -B ^ Xse~^e~ can provide some information on the physics, which occur above the scale ~ lOOGeV. 



2. The B XsC^e" decay has been studied within the standard model and it extensions in 
and ||lO[-Il2| respectively. In the latter works it has been shown that new physics contribution can 
make B — ^ XsC^e" branching ratio two and more times larger than in the SM. These calculations 
have been performed in the next-to-leading order (NLO) of the perturbation theory, which includes 
0{a~^) and 0(1) contribution to i? — Xse^e~ decay. On the other hand, it has been proven in 
that to this order the obtained results may suffer from the uncertainties, connected with those in the 
choice of the low-energy scale ~ and the heavy mass (matching) scale ~ ^t- 

Recently the 0{as) corrections to the B XgC^e^ decay branching have been calculated in 



131 , |I4[. The calculations have been performed in the off-resonance region, corresponding to 0.05 < 
s < 0.25, where s is the dilepton invariant mass, normalized over the b - quark mass. It has been 
found that within the standard model these corrections reduce the above-mentioned uncertainties 
about two times. 

The aim of the present paper is studying the impact of the 0{as) corrections on the B —>■ X^e^e" 
decay rate in the two-Higgs doublet extension of the standard model (2HDM). The investigation is 
carried out for the most general version of the 2IIDM (so-called Model III), where sizable deviations 
from the SM result are possible. The deviations from the SM predictions (called new physics effects) 
occur due to diagrams with the charged Higgs boson mediated loops. It is shown that when the 
contribution of these diagrams is sizable, 0{as) contribution suppresses the B —>■ XsC+e" decay 
width up to 1.5 -7- 3 times, depending on the choice of the dilepton invariant mass and the low - 
energy scale. As a result, in the experimentally allowed region of the parameter space the behavior 
of the B Xse~^e~ branching ratio with the new physics parameters is changed drastically. It is 
also pointed out that the dependence of B{B Xse~^e~) on the parameters of theory can further be 
modified. While due to the 0{as) corrections low - energy scale dependence of the obtained results is 
very small (or even negligible), the matching scale dependence remains large enough. This indicates 
that higher order effects in the perturbation theory can be numerically relevant as well. 



3. It is known that in the SM the decay B —>■ Xge^e" is loop-induced: in the lowest order it 
proceeds via exchange of the up-type quarks and boson in the loops. In the 2IIDM there are 
additional diagrams with W-boson replaced by the charged Higgs boson {H^). The interaction of 
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the charged Higgs boson with quarks may be written in the following form [115 
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nif, are the running top and 
T6[]), and the 



where V is the CKM matrix, g is the weak coupling constant, rrit 
bottom masses (for the relation between the pole and running quark masses see e.g. 
parameters C,t and C,b are the functions of the Higgs doublet vev's, top and bottom masses and the 
couplings of Yukawa interaction of t- and b- quarks with the Higgs doublets fl^. Due to the Higgs 
doublet vacuum phase, C,t and C,b are complex in general. Notice also that while ~ 1 or smaller, 
\C,b\ can be much larger than unity, unless it contradicts with the experimental constraints on the 
B Xg'y branching. 

To avoid flavor changing neutral currents (FCNC) in the Lagrangian, one usually considers the 
simplified versions of the 2HDM: Model I, where only one of the Higgs doublets interacts with 
quarks and the Model II, where one of the Higgs doublets interacts with the up-type quarks and the 
second one does with the down-type quarks. In these versions of the 2HDM, = C,b = — cot jS and 
C,t = — cot j3, = tan/3 respectively (tan/3 is the Higgs vev's ratio). Notice however that due to the 
stringent constraints on the new physics parameters, in the Models I and II the predictions for the 
B Xse~^e~ decay branching coincide (with (10 -r- 15)% accuracy) to those of the standard model. 

In the present paper the most general version of the 2HDM (Model HI), where both of the Higgs 
doublets interact with both the up-type and down-type quarks, is considered (neglecting possible 
FCNC's in the Lagrangian). In this case, due to larger parameter space, B —>■ Xse^e~ branching 
can be up to three times larger than in the SM. It is worth also to mention that the results derived 
for the Model HI may be also considered valid for multi-Higgs doublet models with only one light 
charged Higgs boson. 



4. The B Xse~^e~ decay is studied, using the effective theory with five quarks obtained by 
integrating out the heavy degrees of freedom which are W and Z bosons, t-quark and the charged 
Higgs boson. The effective Hamiltonian for the decay B Xse^e~ can be written as 
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where A| = V^gVtb, = V*^yub, Ci are the coefficients of the Wilson expansion and the full set of 



operators Oi can be found elsewhere [0, As it was mentioned above, current study includes 
only so-called off-resonance region, where 0.05 < s < 0.25. In this case B — >■ Xse~^e~ decay is well 
described by 6 ^ se~^e~ and b —>■ se^e^g partonic transitions. The calculation of these partonic 
transitions includes the following three steps: 

1. The Wilson coefficients Ci at the heavy mass scale, fiw ~ Mw, mt, must be calculated, 
matching the effective and full theories. In the next-to-next-to-leading order (NNLO) the 



matching has to be done at the 0{as) level, i.e. Ci{fiw) = Ci^\f^w) + tts/(4vr)Ci^''(/ivy)- When 
determining matching conditions for the Wilson coefficients one must, generally speaking, take 
into account the difference between the electroweak breaking scale and new physics scale. In 
the 2HDM, Model HI, such a problem occurs when mjj+ ^ nw. Such values of the charged 
Higgs mass are out of the scope of this paper: here one takes mH+ = lOOGeV, 200GeV and 
400GeV. 
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2. The renormalization group equations (RGE) must be used to obtain the Wilson coefficients 
at the low-energy scale fi ~ m^. In the next-to-next-to-leading order this step requires the 
knowledge of the anomalous dimension matrix up to the order a^. 

3. The matrix elements of the operators Oi for the processes b —>■ se+e~, b — > se^e~ g have to be 
calculated. 

Using this procedure, one finds that in the standard model the differential branching ratio of i? — 



XgC'^e decay is given by the following expression [14 



4 (i + ?) \ci;ilM\^ + uRe (cii{'Mc:;i{M))B„ (4) 



where agm is the electromagnetic coupling constant, Bgi is the B — >■ X^e^z/ (experimental) branching 
ratio, z is the ratio of c- and b- quark masses squared and the functions g(z) and k(z) are given in 
131 . The effective new Wilson coefficients are related with the old ones, given e.g. in [||, |13], 0, asQ 
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where the function F/, cuy, cug, cuyg are given in |T^, |^ and the relevant dimension six operators Oi 
are the following: 

O2 = slI^cl CLlybh, = slI^ul ULl^^bi, (6) 

Or = -^m,{fi) SLa'^-'bR F^,, O^ = -^^m,{ii) SLa^'^T^R G%. 

e^ _ _ _ _ 

^9 = TTr^^LYbLe'y^.e, do = -—^SL'y^bie'yf.e 

In the expression for C^llew^ilnew the O(a^) terms have been discarded. Similarly, in (4) for |Cj'^{{^P, 
i=7,9,10, only the terms linear in are retained, when performing numerical calculations. 

0{a^ terms in (6) arise due to the 0{a^ corrections to the Wilson coefficients (the terms pro- 
portional to C'f^^^^ 13); due to the 0{a.^ corrections to < e^e" s\Oi^ > matrix elements and due 

^Our notation for the Wilson coefficients is different from that in rcf. here the superscript denotes 0{aj^) 
and 0(1) contribution to d, whereas always denotes 0{as) contribution. 
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to the singular graphs connected with the gluon bremsstrahlung (6 — > se^e~g) effects |Tj]. The 
contribution of nonsingular bremsstrahlung diagrams and 0{as) corrections to the Oiu matrix 
elements are still unknown. It is expected however for them to be small. 

Both for the SM and the 2HDM the content of the operators Oi is the same, so is the anomalous 
dimension matrix. The new physics effects enter only through the matching conditions. In other 
words, formulae (4) and (5) are valid also in the 2HDM, however for the Wilson coefficients one has 
now 
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The full set of formulae for C^^^ {(j) (so is for C^^{ij) in the chosen operator basis) is given in [|T£ 
Here it is worth to notify only that the contribution of the second term in r.h.s. of (8) dominates 
over the last termQ. For (/x) one has 



Cll^\s) - Cll^\s) (10) 

where C'g^o (s), C'g.'io {^) ^'^d ACq-^-^{s) are given in [0. For Cq^Iq^^ (/i) one can easy deduce that 
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where fj = as(/ivy)/Q;s(/i), (/ivy) are given in C^'' (fiw) is given in |]T9[ and the "magic 

numbers" qf~^^ can be found in |[T3| . 

are unknown yet. For this reason here the 



Unfortunately the matching conditions for Cq iq 



SM values Cg {g (yu) will be used. 

5. The effective Wilson coefficient^ deviate from their SM values differently. Thus Cg^^ devi- 
ates for Cg only by few percents. On the contrary, |C^o | can be about 1.4 times larger than in 
the SMQ (Table 1). 

The largest deviations from the SM interval occur for Cj^"^ . It is known that in the standard model 
extensions Re{C^^^) may have the sign opposite to that in the SM and furthermore unlike the SM the 
dispersive part of C^^^ can be complex (recall that in the SM cf'^'^^^ is real and the imaginary part 

^Moreover, the last term in (8) is about one order smaller than the SM contribution. 

■^Because of the smallness of c'f^'^^^ and condition (7) the qualitative discussion throughout the paper is valid both 
for the old and the new Wilson coefficients, unless the difference is specially notified. 

''One may subsequently expect that omitted here 0(as) corrections to C^^^ will not exceed (fO-f- f5)% and those 
to Cq^^ will be negligible. 
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Table 1: Rw = CiJ^'/Cj? 



mH+ 


lOOGeV 


200GeV 


400GeV 


Rio 


0.36^0.39 


0.23^0.26 


0.12 -^ 0.14 



Table 2: The restrictions on \C^jnew{f^)\fs=o = lOOGeV, coming from the condition B{B 

X,7) = (2.96 ± 0.35) x 10^^ 







fiw = m 




H = mfe/2 


/i = mfe/2 


/i = 2mfe 


/i = mi,/2 


fi = mfe/2 


/i = 2mb 


\(^^ff |2 

1 ^7, new \ 


-.04 ^ 0.13 


.03 ^ .12 


.05 ^ .12 


.05 ^ .13 


.05 ^ .12 


.06 ^ .12 



of C^^Jew arises only due to absorptive parts of the Oi and O2 matrix elements). As for \Cj^^\, it is 
bounded due to the experimental constraints on the B — > branching. In this paper the numerical 
calculations are performed both including and neglecting 0{as) effects, and the derived results are 
compared to each other. When neglecting 0{as) corrections to B{B X^e^e") is it reasonable to 
do the same also for B{B — > ^^7). In this case more conservative bound than (1) should be used: 
one takes therefore [|l|] 1 x 10"^ < B{B X^-f) < 4.2 x 10"^ This gives 0.04 < |Cf ^'^-^p < 0.18 in 
the next-to-leading order. When including 0{as) corrections, the condition (1) puts some constraint 
on \Cy^Jew\^ at the point s = 0. Though this point is out of the consideration, due to the weak 
dependence of C^^^ew 011 s (see formulae of fl^) the aforesaid constraint is essential also for the 



considered range of the dilepton invariant mass. 

Because of discarding 0{a1) terms in the expression for jCy^^^p, there are some problems, 
connected with the bounds on this quantity. As on can see from the Table 2, for ^y/ = m\Y the 
restrictions on jCy^^g^H^Q are highly sensitive to the choice of the low-energy scale. Moreover, for 
fi = mf,/2 the condition (1) allows negative values of \Cj^Jf,J\'^. The unnatural negative values of 
\Cj^iew\^ indicate on the fact that in some regions of the 2HDM parameter space the discarded 0{a1) 
terms are numerically relevant and may not be neglected. Such a situation occurs in particular in the 
case when \cf''^^^^\ ~ 1, however due to the 0{as) corrections \C^^J^J\ ~ \Cj^Jew I even smaller. 
Then it is possible that although T(h ^ 57) ~ \C'i^,iew\^ < 0; the long-distance 0(l/mf) corrections 
(which to the considered order are taken proportional to \C'f'^^^^\'^ [0) drive B Xs^y branching 
to the experimentally allowed interval. 

Of course there is no sense to consider such unphysical possibilities, which arise only due to 
neglecting of higher order corrections to the B —>■ Xgj branching. One can easy avoid to consider 
such unfavorable regions of the parameter space. As one can see from the Table 2, for fiw = the 
reliability of the restrictions on \Cj^^new\^ increases. In particular, the allowed interval of \Cj^J^^\^g^Q is 
weakly sensitive to the choice of the low-energy scale and one may take now 0.05 < \C^lnew\^=a — 0-13 
(this bound turns to be valid also for m//+ = 200GeV and 'mH+ = 400GeV). This result is not sur- 
prising: it is known that for = rrit, 0{as) corrections to the B — >■ Xg'j branching are minimized, 
as compared to the case of fiw = [0 • Thus the above-mentioned problem is avoided, if calcula- 
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tions are performed only for /j^w = rrit. This strategy will be used throughout this paper in the most 
of cases. However, the estimation of possible inaccuracy in the obtained results, connected with the 
uncertainty in the choice of the matching scale, will be done as well. 



6. Let us proceed to the numerical results. During the numerical analysis we use Mw = 80GeV, 
^poie ^ ^^^^ ^ 5)GeV, asiMz) = 0.119 ± 0.002 ||, ml°^' = (4.8 ± 0.15)GeV, mjmi, = 0.29 ± 0.02, 



«em = 1/(130.3 ± 2.3) and Bsi = (10.49 ± 0.46)% (see IjT^ and references therein). The low-energy 
scale fi is varied as /i = mb/2, /i = m^, /i = 2mt,. As it was noted already, the matching scale is 
identified here with the top quark mass. Only when the heavy mass scale dependence is examined, 
the matching scale is chosen as fiw = Mw as well. 

The charged Higgs mass is varied as mH+ = lOOGeV, 200GeV, 400GeV. The restrictions on \^t\ 
and \C,b\ are derived from the requirement for the top and bottom Yukawa couplings to be in the 
perturbativity range in a whole energy interval between the electroweak breaking scale and unification 
scales, and using the experimental constraints coming from the measurements of i? — >■ Xgj branching 



ratio and B — B mixing effects (see the discussion in previous section and IJT^ for more details) 

At the low-energy scale is computed, using its two-loop renormalization group equation [|r 
However, when neglecting 0{as) corrections to B ^ XgC^e' and B Xg'-j branching ratios, one- 
loop result for as(/^) is used. 

In the Wolfenstein parameterization the necessary CKM-factors are given by [^: 

^* - -(l-AV2 + A2(p-zr^)), -^ = X\p-^^)^ (12) 



iKbl 

where A = sin ~ 0.22 and the unitarity triangle parameters p and rj can be obtained from the 
unitarity fits, which yield pi, P3 



p2 + ^2 = 0.423 ± 0.064, 38° < 7 = arctan- < 81° (13) 

P 

The numerical calculations consist of two steps. At first one investigates the renormalization scale 
dependence of the B —>■ Xse~^e~ normalized differential width: 

^^^^ 1 dT{B ^ X^e+e-) 1 dB{B ^ X,e+e-) 



T{B ^ X^e+p) ds Bsi rfs 



The use of R{s) would allow one to compare our results with those of ref.'s [T^, Q. The renormaliza- 
tion scale dependence is examined for some fixed characteristic values of the new physics parameters, 
for the "best fit" values of the CKM parameters {p = 0.19, 77 = 0.37) and for the central values of 
the remaining parameters of theory. 

The next step involves complete investigation for the differential and partially integrated branch- 
ing ratios of B ^ XgC^e^ in the 2HDM and comparison of the obtained results with those in the 
SM. The parameters of theory are varied now in the intervals specified above. 

Let me briefly recall how the situation with the low-energy and matching scale dependence of 
R{s) in the SM looks. On the absence of the 0{as) contribution (then C^jJ^^ in (4) are replaced by 
Cf)^ff ^ i=7,9,10), the /i-dependence of R{s) (in the naive scheme) is ~ 6% or smaller. However, 
such a weak sensitivity of the B X^e+e^ decay width to the choice of the low-energy scale in 
the next-to-leading order is accidental 0. The error from the /i-dependence of R{s) grows up to 
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Table 3: The matching scale dependence of Cio(/u) = C*|q (/i) + as{fi) / (An) C}o (yu) in the SM and in 
the 2HDM for |6| = 1 and mH+ = lOOGeV. 





with 0{as) corrections 


without 0{as) corrections 






Hw = m 


— Mw 


Hw = m 




-4,37 ^ -4.02 


-4.50 -4.10 


-4.78 ^ -4.37 


-4.46 -4.05 




-6.24 ^ -5.64 


-6.17 ^ -5.54 


-6.65 ^ -6.00 


-6.13 ^ -5.50 



13%, when 0{as) corrections to the Wilson coefficients are taken into account ||13|. Only after 0{as) 
corrections to the matrix elements of the operators Oj, j=l,2,7,..,10, are also included, this error is 



reduced down to 6.5% |[T^. As for the error, connected with the uncertainty of the heavy mass scale. 



it is reduced from ~ 15% in the NLO to few percents in the NNLO|. 

The difference between the renormalization scale behavior of R{s) in the SM and 2HDM occurs 
predominantly due to C^^^ . The /i- dependence of C^q^ originates only at 0{as) order and is small 
therefore. There is no essential difference between the nw - dependence of Ciq in the SM and 2HDM 
even for \^t\ ~ 1. As one can see from the Table III, in both of models the nw - error of Ciq is ~ 10% 
in the NLO and ~ (2 - 3)% in the NNLO. 

When investigating the renormalization scale behavior of R{s), the new physics parameters are 
chosen as mH+ = lOOGeV, = 0.5 and a) ^t^b = 6.4, b) ^t^b = 2.8 ± 2.8i. In the case a), the 
dispersive part of Cj^^ is real and has the sign opposite to that in the standard model. In the case 
b), the imaginary part of C^^^ dominates over the real oneQ. 

For the case a), the behavior of the B X^e^e^ normalized differential width as a function of 
the dilepton invariant mass for different choices of the low-energy scale and for fiw = is presented 
in the Fig. 1. As one can see from this figure, in the next-to- leading order the uncertainty of R{s) 
connected with that of the low-energy scale is large enough: it reaches 25% for lower values of s 
and 17% for larger values of s. 0(0;^) corrections reduce the /i-error of R{s) significantly: in the 
naive scheme (where the errors of different terms in (4) are simply summed) it is ~ 10% for lower 
values of s and its sign is flipped, it is only about few percents for larger values of s, and it almost 
disappears for the intermediate values of s. Such a dependence of R{s) on the choice of the low- 
energy scale indicates on the cancellation of /i-errors of the terms in (4). In such cases one usually 
uses so-called Kagan-Neubert method ||2^, 0: in context of the present calculations this approach 
implies the calculation of /i-errors of different blocks of the (effective) Wilson coefficients separately 
and then adding them in quadrature. The result derived in the Kagan-Neubert scheme does not 
differ essentially from that in the naive scheme: now the uncertainty of -R(s), connected with that of 
the low-energy scale, is ~ 12% for lower values of s and ~ 6% for large values of s. This means that 
weak /x-dependence of R{s) in the NNLO is not accidental. 

The investigations in the Kagan-Neubert scheme show also that in the next-to-next-to-leading 
order the main source of the /i-dependence is the term proportional to a{i + |^ \C^^Jew{^A\^ (in the 

^Although this result has been derived in considering 0{as) corrections only to the Wilson coefficients, it is 
easy to check that it remains valid also after taking into account 0{as) corrections to the operators matrix elements. 

^For the cases a) and b) the parameters of theory are within or close to the experimentally allowed values only in 
the NNLO. As it is shown in the section 8, the experimentally allowed intervals of the new physics parameters are 
strongly different in the NLO and the NNLO. 
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Figure 1: Low energy scale dependence of R{s) in the NNLO (solid lines) and NLO (dashed and 
dotted lines) in the 2HDM, case a), for = rrit. The analysis is performed both in the naive and 
the Kagan-Neubert schemes. Unless notified, lower dashed (dotted) and solid lines correspond to 
II — mb/2, middle lines correspond to n = nib and upper lines correspond to /i — 2mi,. 
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This 



NLO there is also sizable contribution from the term proportional to Re {CgjJ*^^{n)C'. 
explains why in the NNLO the low-energy scale dependence of R{s) for larger values of s is especially 
small. 

It is important to stress that in the case a), 0{as) corrections suppress the B X^e+e^ decay 
width about 1.5 -j- 3 times, depending on s and ii. In other words, NNLO corrections to R{s) are 
about (35 65)% of the leading and next-to- leading order terms. It is reasonable therefore to expect 
that 0{al) contribution to the B — > X^e+e" decay width will be nonnegligible as well. This is argued 
also by the investigation of the matching scale dependence of R{s). As one can see from the Fig. 2, 
in the next-to-leading order the uncertainty of R{s) is ~ 15%, when varying the matching scale from 
rrit to Mw- This uncertainty is almost independent on s. In the next-to-next-to-leading order the 
situation is quite different: while for larger values of s, /^ly-error of R{s) decreases down to ~ 10%, 
for lower values of s it increases up to ~ 25%. Such a large error can be reduced only by higher order 
corrections to the B XgC^e' decay width. 

The obtained results for the case b) are presented in the Fig. 3. As one can see from this figure, 
the situation is similar to that of the case a). Due to the 0{as) corrections, the low-energy scale 
dependence of R{s) becomes neghgible, so is the matching scale dependence for larger values of s. 
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R(s)/10-4 




Figure 2: The matching scale dependence of R{s) in the NNLO (sohd hnes) and NLO (dashed and 
dotted hnes) in the case a) for /x = m(,/2. The fat (thin) sohd and dashed (dotted) hnes correspond 
to fj,w = mt {fJ^w = Mw)- 



However for lower values of s, the /ivi/-error of R{s) is of the same order (~ 14%) as without 0{as) 
corrections. 

Thus, although in the 2HDM 0{as) corrections reduce significantly the low-energy scale depen- 
dence of the B — >■ Xse'^e~ decay width, higher order terms in the perturbation theory can still be 
nonnegligible, when the new physics contribution is sizable. The importance of higher order correc- 
tions is manifested by the sensitivity of the obtained results to the choice of the heavy mass scale. 

7. Let us compare now the predictions of the 2HDM for the B XgC^e' branching ratio to 
those of the SM (whole allowed range of the 2HDM parameter space is considered now). The devia- 
tion of B{B Xse~^e~) from the SM results can take place due to the change of the sign of Re{Cj^^) 
(the source I), due to the imaginary part of C^^^ , connected with the Higgs doublet vacuum phase 
(the source II), and due to the deviation of C{q^ from C^l^ (the source III). It is easy deduce from 
the formula (4) that all aforementioned three sources increase the B Xge+e^ branching, once 
Re{C^,ff''') < and \Cllf\ > \. 

The contribution of the source II is expected to be small: it is proportional to[] Irn{Cgl^^f,^)Im{CjjJif,^) 

^During the numerical calculation the terms which make the difference between the B — > X^e+e^ decay and the 
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Figure 3: The same as in Fig.'s 1,2 (respectively the left- and right-hand-side graphs) but for the 
case b) . 



and one has \Im[Cf^^ 
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< \Re{Cff 



9 •new I 



The effects connected with the source 111 are strongly cor- 
related with the sign o^Re{C^^^). In the SM-like case the terms of (4), proportional to ICg^^g^p and 
R^iCgfnewCj^Jew) 1 havc oppositc sigus SO that their contribution is partially canceled. Consequently, 
the partial weight of the the term proportional to the [C^o^g^p and hence the contribution of the 
source III is significant. Respectively, when i?e(Cy^^) has the sign, opposite to that in the standard 
model, the effects of the source III are minimized. It is straightforward to deduce from current 
discussion that the sources 1 and 111 do not interfere in fact. 

The maximum value of the B XsC^e" differential branching ratio as a function of s is presented 
in the Fig. 4 for particular cases, when only one of the aforementioned sources is actual. As it was 
expected, the contribution of the source II is not large: if only this source is actual, dB{B — > 
Xse'^e~)/ds deviates from its SM maximum value at most 1.4 times. The source III can make 
dB(B Xse~^e~)/ds 1.6 1.9 times larger than in the SM (respectively for s varying from 0.05 to 
0.25). The deviations from the SM results, connected with the source 111, are most perceptible when 
\^b\ 1 \^t\ ~ 1- The effects of the source 111 are rapidly minimized with the increasing of the 
charged Higgs mass: for m//+ = 400GeV they are already of the same order as those connected with 
the source II. 



CP-conjugated decay are dropped. 
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ds 




Q g I I I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 

s 

Figure 4: Maximum value of dB{B — * Xse~^e~)/ds in the SM (dashed hne) and 2HDM (solid 
hues). The following regions of the 2HDM parameter space are considered: /m(^(^{,) = 0, 

|Croie«,P < l-15|Ci^{e«; P (Imc 1, Only the source I is actual); sign{Re{C^^J^^J) = sign{Re{Cl^J^^ )), 
|Citne«,l' < l-15|C^foTei'l' (^^^^ 2, ouly the source II is actual); Im{Ub) = 0, stgn{Re{C^jJ,J) = 
sign{Re{Cj^^^^), mH+ — lOOGeV, 400GeV (hnes 3, 3' respectively, only the source III is actual). 

The most sizable deviations from the SM predictions occur due to the source I. Due to this source, 
maximum value of dB{B — > Xse'^e'~)/ds can be up to 2.6 times larger than in the standard model. 
The effects connected with the source I are almost insensitive to the charged Higgs mass. Large 
deviations from the SM result are possible even for ~ ITeV (until Re[Cj^^) becomes positive, 
when |(^{,| 3> 1). However the consideration of such large values of the charged Higgs mass is out of 
the scope of the present paper, otherwise one should take into account the difference between the 
electroweak breaking scale and the charged Higgs mass scale. 

It is worth to make here some digression on the situation in the Models I and II. If one applies 
the restrictions on the parameters and notified in the section 6, only the source III will be 
relevant in these versions of the 2HDM. But the effects of this source are very small, because one has 
iCtp < 1/4 for mjy+ = lOOGeV in the Model I and mH+ > 400GeV in the Model II. As a result, the 
predictions of these models for B — > Xse^e~ branching are close to those of the SM. 

Let us return back to the Model III. The results derived for the most general case (when all of 
three sources are actual) are presented in the Fig. 5. As one can see from this figure, while the 
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ds 



^/10-5 
ds 




0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 



s s 

Figure 5: Minimum and maximum values of dB{B — > Xse^e~)/ds in the SM (dashed hnes) and 
2HDM (sohd hnes) a) in the NNLO, b) in the NLO. 



minimum value of dB{B — > XgC'^e') / ds coincides (within 10% accuracy) with that in the SM, the 
maximum value is 2.5 3 times larger than in the standard model. Note that 0(0;^) corrections do 
not enhance or suppress the deviation from the SM predictions. This is because in both the NLO 
and the NNLO the experimental constraints on the B — > Xg^ branching fix |Cy'^''^| approximately in 
the same range (at least for the choice of the heavy mass scale jiw — n^t)- 
For the partially integrated branching ratio 

^ ' Jo.05 ds ^ ' 

one gets: 

SM, NLO: AS(S ^ X.e+e") = (1.3 ^ 2.1) x 10"^ 
2HDM, NLO: AS(S ^ X^e+e") = (1.3 ^ 5.6) x 10"^ 
SM, NLO: A5(B ^ X,e+e-) = (1.2 ^ 1.9) x 10"^ 
2HDM, NLO: A5(5 ^ X.e+e") = (1.1 ^ 5.3) x 10"^ 

As it follows from these results, in the 2HDM partially integrated branching ratio of i? — > XgC^e" 
decay can be up to 2.8 times larger than in the standard model. Again, the deviations from the SM 
result are not affected by 0{as) corrections. 

Thus, if both B — > and B — > Xse^e~ branching ratios are computed with the same accuracy, 
then in any order of the perturbation theory (unless in some numerically relevant order C^^^ becomes 
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Figure 6: Minimum and maximum values of the partially integrated branching ratio of -B ^ XgC+e" 
decay as functions of \itib\ for mH+ = 400GeV in the NLO (dashed lines) and NNLO (fat solid 
lines). Straight thin solid and dotted hnes represent the SM interval respectively in the NNLO and 
the NLO. 



highly sensitive to s) one will derive that in the 2HDM, maximum value of the B — > XgC'^e branching 
ratio is 2.5 -^ 3 times larger than in the standard model. 

8. When considering whole allowed range of the 2HDM parameter space, one can find out the 
difference between the NLO and NNLO predictions for the B — > Xse'^e~ branching, examining 
the dependence of this variable on the new physics parameters. Once in the 2HDM 0(0;^) terms are 
important primarily for C^^^ew^ is reasonable to take mH+ = 400GeV, to minimize the contribution 
of the source III, and investigate the dependence of the B — > XgC^e' branching ratio on the product 
\itih\ (see formula (8) and the comment beneath). The minimum and maximum values of the partially 
integrated B — > Xge^e' branching ratio as functions of \itib\ are presented in the Fig. 6. One can see 
that the dependence of AS on \^t^b\ is quite different in the NLO and the NNLO. Thus, in the next-to- 
leading order maximum value of AB{B X^e^e") is two and more times larger than in the standard 
model, when taking 8 < \it^b\ < 15. In the next-to-next-to-leading this occurs for 17 < \^t^b\ < 30. 
In other words, after including 0{as) corrections, the dependence of AB{B — > X^e+e") on \^tib\ is 
changed drastically. 
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Generally speaking, the behavior of AB{B X^e+e") with \C,t^b\ can further be modified by 
higher order effects in the perturbation theory. On the other hand, it has been shown that in the 
next-to- next-to-leading order, when choosing the matching scale as fiw = ^t, the low-energy scale 
dependence of the B X^e^e" decay width is weak. In order to preserve the weak sensitivity of 
R{s) to n, higher order corrections to the B Xge^e" branching must somehow cancel each other. 
This allows one to expect that higher order effects in the perturbation theory will not modify the 
obtained result so drastically, as it is the case in the Fig. 6. In other words, one may suppose that 
the dependence of AB{B — > XsC'^e') on derived in the NNLO, is close enough to the proper 

one. However, it is necessary to emphasize that such a suggestion may be done only for specific 
choices of the heavy mass scale (like above). More generally (for instance for fiw = Mw) the NNLO 
results^ for the B — > Xse^e~ branching can be highly unreliable and be largely modified by higher 
order effects in the perturbation theory. 

Summarizing the discussion of this section one may conclude that including the 0{as) effects is an 
important step on the way of establishing the proper relations between the B — > X<.e"'"e~ branching 
and the new physics parameters. 

9. Thus, 0{as) corrections to B ^ X^e'^e' decay have been examined in the two-Higgs dou- 
blet extension of the standard model. The investigations have been performed for the most general 
version of the 2HDM (Model III) in the so-called off-resonance region of the dilepton invariant mass 
(0.05 < s < 0.25). 

It has been shown that in the case when the new physics effects are sizable, 0(0;^) corrections 
(for the fixed values of the 2HDM parameters) suppress the B — > Xse~^e~ decay width 1.5 ^3 times. 
It is natural to suppose that O(a^) corrections to the B — > Xse~^e~ branching will be numerically 
relevant as well. 

Last suggestion is confirmed by the fact that the obtained results are still sensitive to the choice 
of the heavy mass scale, when the new physics effects are sizable. Even after 0{as) corrections are 
included, the uncertainty in the B Xse^e~ decay width connected with the choice of the matching 
scale reaches 25%. 

On the other hand, 0{as) corrections reduce significantly the low - energy scale dependence of 
the B — >■ Xse^e~ branching. The /i - error of the obtained results in the NNLO is ~ 10% or smaller 
(compared to ~ 25% in the NLO). This means that in the next-to- next-to leading order the reliability 
of the predictions for the B XgC^e^ branching ratio increases. 

When using the experimental constraints on the B Xg'y branching and calculating B X^j 
and B —>■ X^e+e" decays with the same accuracy, one will probably get in all orders of the pertur- 
bation theory that in the 2HDM B —>■ X^e^e" branching ratio can be about 2.5 -r- 3 times larger 
than in the standard model. However only after taking into account 0{as) corrections and probably 
those of higher orders, one is able to derive the proper relations between the new physics parameters 
and the B — X^e^e" branching ratio. 

Author is grateful to H. M. Asatrian and H.H. Asatryan for stimulating discussions. This work 
has been partially supported by the SCOPES program. 

^When taking = Mw, one can examine the dependence of the B — > X^e'^e^ branching on the new physics 
parameters, requiring that < | C'f{/eu) (m) ll=o — 0-13- 
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